Transition from single-file to two-dimensional diffusion of interacting particles in a 

quasi-one-dimensional channel 
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Diffusive properties of a monodisperse system of interacting particles confined to a quast-one- 
dimensional (QID) cliannel are studied using molecular dynamics (MD) simulations. We calculate 
numerically the mean-squared displacement (MSD) and investigate the influence of the width of 
the channel (or the strength of the confinement potential) on diffusion in finite-size channels of 
different shapes (i.e., straight and circular). The transition from single-file diffusion (SFD) to the 
two-dimensional diffusion regime is investigated. This transition (regarding the calculation of the 
scaling exponent (a) of the MSD {Ax^{t)) oc t") as a function of the width of the channel, is shown 
to change depending on the channel's confinement profile. In particular the transition can be either 
smooth (i.e., for a parabolic confinement potential) or rather sharp/stepwise (i.e., for a hard-wall 
potential), as distinct from infinite channels where this transition is abrupt. This result can be 
explained by qualitatively different distributions of the particle density for the different confinement 
potentials. 



PACS numbers: 05.40.-a, 66.10.C-, 82.70.Dd, 83.10.Rs 



I. INTRODUCTION 

There is a considerable theoretical and practical inter- 
est in the dynamics of systems of interacting particles 
in confined geometries Single- file diffusion (SFD) 

refers to a one-dimensional (ID) process where the mo- 
tion of particles in a narrow channel (e.g., quasi-lD sys- 
tems) is limited such that particles are not able to cross 
each other. As a consequence, the system diffuses as 
a whole resulting in anomalous diffusion. The mecha- 
nism of SFD was first proposed by Hodgkin and Keynes 
in order to study the passage of molecules through 
narrow pores. Since the order of the particles is con- 
served over time, this results in unusual dynamics of 
the system 0, H], different from what is predicted from 
diffusion governed by Pick's law. The main character- 
istic of the SFD phenomena is that, in the long-time 
limit, the MSD (mean-square displacement, defined as 

{Ax^t)) = {EZii^/N)[x^it + At) ~ x,it)]^)At) scales 
with time as 



0.5 



(1) 



This relation was first obtained analytically in the pi- 
oneering work of Harris sj. Recent advances in nan- 
otcchnology have stimulated a growing interest in SFD, 
in particular, in the study of transport in nanopores 
Ion channels of biological membranes and carbon nan- 
otubes [1] are examples of such nanopores. The macro- 
scopic flux of particles through such nanopores is of great 
importance for many practical applications, e.g., parti- 
cle transport across membranes is a crucial intermediate 
step in almost all biological and chemical engineering pro- 
cesses. SFD was observed in experiments on diffusion of 
molecules in zeolite molecular sieves 0. Zeolites with 
unconnected parallel channels may serve as a good real- 



ization of the theoretically investigated one-dimensional 
systems. SFD is also related to growth phenomena |10j. 

The theoretical background of SFD was developed in 
early studies on transport phenomena in ID channels [TTI - 
[Tsj . It is also interesting to learn how the size of the sys- 
tem will influence the diffusive properties of the system. 
SFD in finite size systems has been the focus of increas- 
ing attention since there are few exact theoretical results 

to date EMI, 

which showed the existence of different 
regimes of diffusion. 

Colloidal systems, complex plasmas and vortex mat- 
ter in type-II superconductors are examples of systems 
where SFD may occur. The use of colloidal particles is 
technically interesting since it allows real time and spa- 
tial direct observation of their position, which is a great 
advantage as compared to atoms or molecules, as shown 
recently in, e.g., the experimental study of defect induced 
melting |17| . One typically uses micro-meter size col- 
loidal particles in narrow channels, as shown in (isl. [l9|. 
The paramagnetic colloidal spheres of 3.6 fim were con- 
fined in circular trenches fabricated by photolithography 
and their trajectories were followed over long periods of 
time. Several other studies have focused on the diffusive 
properties of complex plasmas. A complex plasma con- 
sists of micrometer-sized ("dust") particles immersed in 
a gaseous plasma background. Dust particles typically 
acquire a negative charge of several thousand elementary 
charges, and thus they interact with each other through 
their strong electrostatic repulsion [20j . 

Systems of particles moving in space of reduced dimen- 
sionality or submitted to an external confinement poten- 
tial exhibit different behavior from their free-of-border 
counterparts [2l[ . The combined effect of interaction be- 
tween the particles and the confinement potential plays a 
crucial role in their physical and chemical properties j22| . 
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In Ref. [l^l , it was found that SFD depends on the inter- 
particle interaction and can even be suppressed if the in- 
teraction is sufficiently strong, resulting in a slower sub- 
diffusive behavior, where (Ax^(t)) cx t", with a < 0.5. 

In this paper, we will investigate the effects of confine- 
ment potential on the diffusive properties of a QID sys- 
tem of interacting particles. In the limiting case of very 
narrow (wide) channels, particle diffusion can be referred 
to SFD (2D regime) characterized by a subdiffusive (nor- 
mal diffusive) long-time regime where the mean-squared 
displacement (MSD) (Ax^{t)) oc t°-^ (oc t^ °). Recah that 
the MSD of a tagged hard-sphere particle in a one dimen- 
sional infinite system is characterized by two limiting dif- 
fusion behaviors: for time scales shorter than a certain 
crossover time = where D is the diffusion coef- 

ficient and p is the particle concentration, {Ax'^{t)) cx t^-° 
which is referred to as the normal diffusion regime [23 |. 
For times larger than Tc, the system exhibits a subdiffu- 
sive behavior, with the MSD {Ax'^{t)) cx t*^'^, which char- 
acterizes the single-file diffusion regime. Between these 
two regimes, there is a transient regime exhibiting a non- 
trivial functional form. 

However, in case of a finite system of diffusing particles 
(e.g., a circular chain or a straight chain in the presence of 
periodic boundary conditions), the SFD regime (i.e., with 
(Ax^(<)) cx t'^-^) does not hold for t oo, unlike in an 
infinite system. Instead, for sufficiently long times, the 
SFD regime turns to the regime of collective diffusion, 
i.e., when the whole system diffuses as a single "parti- 
cle" with a renormalized mass. This diffusive behavior 
has been revealed in experiments [l^, and theoreti- 
cal studies [2^, [26l - [29| . This collective diffusion regime is 
similar to the initial short-time diffusion regime and it is 
characterized by either (Ax'^(t)) cx t^ '^, for overdamped 
particles (see, e.g., [Hill, |2|]) or by {Ax'^{t)} cx t^-° 
(followed by the MSD oc t^ °), for underdamped systems 
[27l [29| . Correspondingly, the time interval where the 
SFD regime is observed becomes finite in finite size sys- 
tems. It depends on the lenght of the chain of diffusing 
particles: the longer the chain the longer the SFD time 
interval. Therefore, in order to observe a clear power- 
law behavior (i.e., {Ax'^{t)) (x t") one should consider 
sufficiently large systems. 

Here we focus on this intermediate diffusion regime 
and we show that it can be characterized by {Ax^{t)) 
oc t", where 0.5 < a < 1.0, depending on the width (or 
the strenght of the confinement potential) of the chan- 
nel. We analyze the MSD for two different channel ge- 
ometries: (i) a linear channel, and (ii) a circular channel. 
These two systems correspond to different experimental 
realizations of diffusion of charged particles in narrow 
channels [3 The latter one (i.e., a circular chan- 
nel) has obvious advantages: (i) it allows a long-time 
observation of diffusion using a relatively short circuit, 
and (ii) it provides constant average particle density and 
absence of density gradients (which occur in, e.g., a lin- 
ear channel due to the entry/exit of particles in/from 
the channel). Thus circular narrow channels were used 



in diffusion experiments with colloids [18| and metallic 
charged particles (balls) Furthermore, using differ- 
ent systems allows us to demonstrate that the results 
obtained in our study are generic and do not depend on 
the specific experimental set-up. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model and numerical approach. In Sec. HI 
diffusion in a system of interacting particles, confined to 
a straight hard-wall or parabolic channel, is studied as 
a function of the channel width or confinement strength. 
In Sec. IV, we discuss the possibility of experimental ob- 
servation of the studied crossover from the SFD to 2D 
diffusive regime. For that purpose, we analyze diffusion 
in a realistic experimental set-up, i.e., diffusion of mas- 
sive metallic balls embedded in a circular channel with 
parabolic confinement whose strength can be controlled 
by an applied electric field. The long-time limit is ana- 
lyzed in Sec. V using a discrete site model. Finally, the 
conclusions are presented in Sec. VI. 

II. MODEL SYSTEM AND NUMERICAL 
APPROACH 

Our model system consists of N identical charged 
particles interacting through a repulsive pair potential 
Vi„t(rij). In this study, we use a screened Coulomb po- 
tential (Yukawa potential), Vint exp(— r/A£))/r. In 
the transverse direction, the motion of the particles is 
restricted either by a hard-wall or by a parabolic con- 
finement potential. Thus the total potential energy of 
the system can be written as: 

N N 
i—1 i>j — l 

The first term in the right-hand side (r.h.s.) of Eq. ([2]) 
represents the confinement potential, where Vcir'i) is 
given by: 

Vc{r.) = l^ f'l^i-^-il (3) 
" [00 for \y,\ > ^ ' 

for the hard- wall confinement, 

V{n) = ^mu^yl (4) 

for parabolic one-dimensional potential (in the y- 
direction), and by 

Vein) ^ Piro - (5) 

for parabolic circular confinement. Here is the width 
of the channel (for the hard- wall potential) , m is the mass 
of the particles, ujq is the strenght of the parabolic ID 
confining potential, rg is the coordinate of the minimum 
of the potential energy and is the displacement of the 
ith particle from tq (for the parabolic circular potential) . 
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Note that in case of a circular channel, tq = rch, where 
r^h is the radius of the channel. 

The second term in the r.h.s. of Eq. ([2]) represents 
the interaction potential between the particles. For the 
screened Couloumb potential, 



Vintinj) 



q2 g-ln-fjl/Ao 



(6) 



where q is the charge of each particle, e is the dieletric 
constant of the medium, — \ri~fj \ is the distance be- 
tween ith and jth particles, and is the Debye screen- 
ing length. Substituting ^ into Eq. ([2]), we obtain the 
potential energy of the system Hy '- 



N 



(7) 



i>]=i 



In order to reveal important parameters which character- 
ize the system, we rewrite the energy Hy in a dimension- 
less (Hy) form by making use of the following variable 
transformations: Hy — (<7^/eao)-ffy , r = r'ag, where uq 
is the mean inter-particle distance. The energy of the 
system then becomes 



N 



N 



(8) 



4=1 



l>j = l 



where k — uq/ Xjj is the screening parameter of the in- 
teraction potential. In our simulations in Sec. Ill, we 
use a typical value of k = 1.0 for colloidal systems and 
Xd = lO-^m. 

The hard-wall confinement potential is written as 



for ly.'l < i?^/2 
oo for \y'^\ > 



(9) 



In Eq. dTTI) . ri, Di and are the position, the self- 
diffusion coefficient (measured in m^/s) and the mass (in 
kg) of the ith particle, respectively, t is the time (in sec- 
onds), fcs is the Boltzmann constant, and T is the abso- 
lute temperature of the system. Finally, is a randomly 
fluctuating force, which obeys the following conditions: 
{Ft) = and {F^{t)F^{t')) = 2rikBT6,i'S{t - t'), where 
rj is the friction coefficient. Eq. can be written in 

dimensionless form as follows: 



dr', 
If 



(12) 



where we use the following transformation Vint — 
(q^/eao)T^„j, D[ — Di/a^, and introduced a coupling pa- 
rameter r, which is the ratio of the average potential 
energy to the average kinetic energy, F = (V) / {K) , such 
that F = /ksTeao. The time t' is expressed in seconds 
and distances are expressed in units of the interparticle 
distance oq. In what follows, we will abandon the prime 
(') notation. We have used a first order finite difference 
method (Euler method) to integrate Eq. ([T^ numerically. 
In the case of a straight channel, periodic boundary con- 
ditions (PBC) were applied in the x-direction while in 
the y-direction the system is confined either by a hard- 
wall or by a parabolic potential. Also, we use a timestep 
At = 0.001 and the coupling parameter is set to F = 10. 
For a circular channel, we use polar coordinates (r, 0) and 
model a 2D narrow channel of radius rch with parabolic 
potential-energy profile across the channel, i.e., in the 
r-direction. 



III. ID VERSUS 2D DIFFUSION IN A 
STRAIGHT CHANNEL 



where is scaled by the inter-particle distance ag. We 
also introduce a dimensionless parameter 



X 



m{uJoaoy 
2kBT 



(10) 



which is a measure of the strenght of the parabolic ID 
confinement potential. 

For colloidal particles moving in a nonmagnetic liq- 
uid, their motion is overdamped and thus the stochastic 
Langevin equations of motion can be reduced to those 
for Brownian particles (sij : 



dt 



kBT\ 



-y^Vc{n) + F^^{t) 



(11) 



Note, however, that in Sec. IV we will deal with massive 
metallic balls and therefore we will keep the inertial term 
in the Langevin equations of motion. 



A. Mean-square displacement (MSD) calculations 

In order to characterize the diffusion of the system, we 
calculate the MSD as follows: 



N 



i=i 



At 



(13) 



where N is the total number of particles and (...)At rep- 
resents a time average over the time interval At. Note 
that in the general case (e.g., for small circular channels 
with the number of particles = 20 — see Sec. IV) 
the calculated MSD was averaged over time and over the 
number of ensembles [s^. However, we found that for 
large N (i.e., several hundred) the calculated MSD for 
various ensemble realizations coincide (with a maximum 
deviation within the thickness of the line representing the 
MSD). 

To keep the inter-particle distance approximately equal 
to unity, we defined the total number of particles N for 
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a ID and QID system as 



N = 



L 



Rw < 1, 



(14) 



where L is the size of the simulation box (in dimensionless 
units) in the x-direction. In our simulations for a straight 
channel geometry, we typically used N = 400 — 900 par- 
ticles. We study the system for two different types of 
confinement potential: (i) a parabolic ID potential in 
the y-direction, which can be tuned by the confinement 
strength x and (ii) a hard-wall potential, where particles 
are confined by two parallel walls separated by a distance 
. The results of calculations of the MSD as a function 
of time for different values of the confinement strength x 
[Eq. (|T0)) ] and the width of the channel are presented 
in Figma)-(c) and FigI51Ja)-(c), respectively. 

Initially, in both cases (i.e., a parabolic and a hard- 
wall confinement potential), the system exhibits a short- 
time normal diffusion behavior, where {Ax'^{t)) oc 
This is the typical initial "free-particle" diffusion regime. 
After this initial regime, there is an intermediate sub- 
diffusive regime (ITR). As discussed in Ref^jSGj, the 
ITR shows an apparent power-law behavior j37| . where 
0.5 < a < 1.0, and it was also found previously in dif- 
ferent diffusion models [38, 39]. In the ITR, we found 
a SFD regime for either a channel with strong parabolic 
confinement [x = 3.5 (Fig. [ija))] or a narrow hard-wall 
channel 0.20 (Fig. d^a))]. This is due to the fact 

that for large (small) values of x (R-w), the confinement 
prevents particles from passing each other. The results 
for a in the ITR are shown as a function of x a-nd R^, in 
Fig. [Hd) and Fig. [DJd), respectively. As can be seen in 
Fig. Hid) [Fig.lDJd)], a increases with decreasing x [with 
increasing R^] and thus the SFD condition turns out to 
be broken. The values of a presented in these figures cor- 
respond to the minimum of the effective time dependent 
exponent a{t). Following Ref. [i^, a{t) is calculated us- 
ing the "double logarithmic time derivative" 



ait) 



d log(Aa;^(t)) 
d log t 



(15) 



and the results are shown in Fig. [3l 

The different diffusive regimes, i.e. normal diffusion 
regime (a = 1.0) and SFD (a = 0.5), were also found re- 
cently in finite-size systems [13, [H| although the transi- 
tion from SFD to normal diffusion was not analyzed. The 
a-dependence on both the confinement parameters (i.e., 
a{x) and Q;(i?u,)) presents a different qualitative behav- 
ior, namely, the SFD regime is reached after a smoother 
crossover in the parabolic confinement case as compared 
to the hardwall case. A similar smoother crossover is 
also found in the case of a circular channel with parabolic 
confinement in the radial direction. A more detailed dis- 
cussion on these two different types of the behavior of a 
will be provided in Sec. IIVI 



10"' 


- (a) 


.1.0 . , 

/ 


10"" 


















J-J 

< 10-'^ 


^0.5 ^^^^^^ 




V 






10" 




r = 10|x = 3.5|: 



10' 10' 10° 10' 10' 10' 10' 

t(s) 



10' 

_ 10'' 
E 

A 10 



V 

10" 



r (b) 








0.62 ^,1^ 










r=10|x = 2.5|- 



10"' 10"' 10° 10' 10' 10^ 10' 

t(s) 



10"' 



10"' 



10"' 



X10"^ 



10"' 



r (c) 








0.78 










r = 10|x = l.5|i 



10"' 10"' 10° 10' 10' 10' 10' 
t(s) 




FIG. 1: (Color online) (a)-(c) Log-log plot of the mean-square 
displacement (MSD) {Aa;^(t)) as a function of time for differ- 
ent values of x- Different diffusion regimes can be distin- 
guished: normal diffusion regime (q = 1.0) and intermediate 
subdiffusive regime (ITR, a < 1.0). Note that for the case of 
X = 1.5, there is a normal diffusion regime (i.e. a = 1.0) after 
the ITR. The dashed and solid lines in (a)-(c) are a guide to 
the eye. Panel (d) shows the dependence of the slope (a) of 
the MSD curves (in the ITR, characterized by an apparent 
power-law; {Ax^{t)) oc t"") on the confinement strength x- 
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FIG. 2: (Color online) (a)-(c) Log-log plot of the mean-square 
displacement (MSD) {Ax^{t)) as a function of time for differ- 
ent values of Rw Different diffusion regimes can be distin- 
guished: normal diffusion regime (a = 1.0) and intermediate 
subdiffusive regime (ITR, a < 1.0). Note that for the case of 
Rn, = 0.60, there is a normal diffusion regime (i.e. a = 1.0) 
after the ITR. The dashed and solid lines in (a)-(c) are a guide 
to the eye. Panel (d) shows the dependence of the slope (a) 
of the MSD curves (in the ITR, characterized by an apparent 
power-law; {Ax'^{t)) oc f°) on the confinement parameter R^. 
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FIG. 3: (Color online) (a)-(b) Exponent a as a function of 
time, calculated from Eq. (|15p for different values of the con- 
finement parameters x and R^, respectively. 



B. "Long-time" behavior of the MSD curves and 
crossing events C{t) 

For small values of the parabolic confinenient (e.g., 
X = 1.5), the MSD curves present three different dif- 
fusive regimes: (i) a short-time normal diffusion regime, 
where MSD {Ax'^{t)) oc t^-°; (ii) a subdiffusive regime 
with {Ax^{t)) (X t", where 0.5 < a < 1.0 and (iii) 
a "long-time" diffusion regime, which is characterized 
by {Ax^{t)) oc t^ °. Note that the "long-time" term 
used here is not to be confused with the long-time used 
for infinite systems, as discussed in the Introduction. 
However, for large values of the parabolic confinement 
(e.g., X — 3.5), we observe only two distinct diffusive 
regimes, namely: (i) a short-time normal diffusion regime 
((Ax2(<)) oc ti O) and (ii) a SFD regime (i.e., {Ax'^it)) 

One question that arises naturally is whether this nor- 
mal diffusion regime (i.e., {Ax'^{t)) oc ° for "long- 
times" ) is an effect of the colletive motion of the system 
(center-of-mass motion) or an effect of the single-particle 
jumping process, since the confinement potential x = 1-5 
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allows particles bypass. In order to answer this question, 
we calculate the number of crossing events C(t) as a func- 
tion of time and results are shown in Fig.lH^a). We found 
that for small values of the confinement potential (e.g., 
X = 1.5) the number of crossing events grows linearly in 
time, i.e., C(t) oc ujct, where Uc is the rate of crossing 
events. On the other hand, a strong confinement poten- 
tial (e.g., X = 3.5) prevents particles from bypassing, and 
thus C(t) = during the whole simulation time. 

Therefore, the "long-time" normal diffusive behavior 
(i.e., {Ax'^{t)) cx t^''^ for "long-times") found in our sim- 
ulations for the case where the SF (single- file) condition is 
broken (e.g., x = 1-5) is not due to a collective (center-of- 
mass) diffusion. Instead, this normal diffusive behaviour 
is due to a single-particle jumping process, which hap- 
pens with a constant rate Wc > for the case of small 
values of the confinement (x — 1.5) and cjc = (for 
X = 3.5). The same analysis was done for the case of 
the hard-wall confinement potential, and the results are 
found to be the same as for the parabolic confinement. 

Nevertheless, we point out that the collective diffusion 
does indeed exist, but our results from simulations do not 
allow us to observe this collective (center-of-mass) diffu- 
sion regime because of the large size of our chain of par- 
ticles {N = 400 - 900). Simulations with TV = 80 - 100, 
and excluding the possibility of mutual bypass (strong 
confinement potential), allowed us to observe that the 
(Aa;^(t)) cx ° regime is recovered in the "long-time" 
limit. In Sec. V, we will further discuss the long-time 
limit using a model of discrete sites. 

As we demonstrated above, the transition from pure 
ID diffusion (SFD) characterized by a = 0.5 to a 
quasi-lD behavior (with a > 0.5) could be either more 
"smooth" (as in Fig. 1(d), for a parabolic confinement) 
or more "abrupt" (as in Fig. 2(d), for a hard- wall con- 
finement). One can intuitively expect that this difference 
in behavior can manifest itself also in the crossing events 
rate uJc, i.e., that ujc as a function of x (or Rw) should 
display a clear signature of either "smooth" or "abrupt" 
behavior. 

However, the link between the two quantities, i.e., 
the exponent, a{x/Rw), and the crossing events rate, 
Wcix/Rw) is not that straightforward. To understand 
this, let us refer to the long-time limit (which will be ad- 
dressed in detail within the discrete-site model in Sec. V). 
As we show, in the long-time limit the exponent a is de- 
fined by one of the two conditions: Uc — Q (then a — 0.5) 
or > (then a — 1) and it does not depend on the 
specific value of uj^ provided it is nonzero. Therefore, in 
the long-time limit the transition between ID to 2D be- 
havior is not sensitive to the particular behavior of the 
function ujc{x/ Rw)- 

Although for "intermediate" times (considered in this 
section) the condition Wc = or cjc > is not critical, 
nevertheless, very small change in the crossing events rate 
uJc{x/Rw) strongly influences the behavior of the expo- 
nent a{x/Rw)- This is illustrated in Figs. 4(b, c). In 
Fig. 4(b), the function Ludx) gradually decreases from 
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FIG. 4: (Color online) (a) Number of crossing events C{t) as 
a function of time for A'' = 400 particles, for x ~ 1-5 (black 
open circles) and x ~ 3.5 (green open diamonds). The solid 
red line is a linear fit to C{t). Panels (b) and (c) show the 
rate of the crossing events cjc as a function of the confinement 
potential parameters (x and Rw)- The insets in the panels (b) 
and (c) show the derivatives, dwc{x)/dx and duuc{Rw) / dR^^ , 
correspondingly. 



1.45 to for X varying in a broad interval from 1.5 to 
3 (note that the segment of Wc(x) for 2.5 < x < 3 
is nonzero which can be seen in the inset of Fig. 4(b) 
showing the derivative dujc{x)/dx)- Correspondingly, the 
transition from a — 0.5 to a « 0.8 in that interval of x 
is "smooth" (see Fig. 1(d)). On the other hand, the 
function ujc{Rw) shown in Fig. 4(c) mainly changes (note 
the change of the slope dwc{Rw) / dRw shown in the inset 
of Fig. 4(c)) in a narrow interval 0.5 < Rw < 0.6. Re- 
spectively, the transition for the function a{Ria) occurs 
in the narrow interval 0.5 < R^j < 0.6 and thus is (more) 
"abrupt" . 



C. Distribution of particles along the t/-direction 

For the ideal ID case, particles are located on a straight 
line. Increasing the width Ryj of the confining channel 
will lead to a zig-zag transition [4l[. This zig-zag 
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FIG. 5: For the hard- wall confinement case, we show typical 
trajectories of particles (i.e. 10^ MD simulation steps) con- 
fined by the channel of width (a) R^^ ~ 0.20, (b) = 0.60 
and (c) R^ = 0.80. 



configuration can be seen as a distorted triangular con- 
figuration in this transition zone. Further increase of 
brings the system into the 2D regime, where the normal 
diffusion behavior is recovered (see Fig. [5]). 

For the parabolic ID confinement, we can see 
[Fig. [6ja)] that the distribution of particles P{y) along 
the channel is symmetric along the axis y = 0. Also, 
for large values of % (e.g., x = 3.5) particles are con- 
fined in the y-direction and thus can move only in the 
x-direction, forming a single-chain structure. As the con- 
finement decreases (x — >■ 0), the distribution of parti- 
cles P{y) broadens resulting in the crossover from the 
SFD regime (% = 3.5) to the 2D normal diffusion regime 
(x = 0.5). Note that for small values of x (e.g., X = 0.5), 
the system forms a two-chain structure (represented by 
two small peaks of P{y) in Fig.jBJa)), thus allowing par- 
ticles to pass each other. 



IV. DIFFUSION IN A CIRCULAR CHANNEL 

In the previous section, we analyzed the transition 
(crossover) from the SFD regime to 2D diffusion in nar- 
row channels of increasing width. The analysis was per- 
formed for a straight channel with either hard-wall or 
parabolic confinement potential. However, in terms of 
possible experimental verification of the studied effect, 
one faces an obvious limitation of this model: although 
easy in simulation, it is hard to experimentally fulfill the 
periodic boundary conditions at the ends of an open chan- 
nel. Therefore, in order to avoid this difficulty, in SFD 
experiments [H, HI] circular channels were used. 

In this section, we investigate the transition (crossover) 
from SFD to 2D-diffusion in a system of interacting par- 
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FIG. 6: (Color online) Probability distribution of the particle 
density P{y) along the j/-direction are shown for (a) different 
values of x (parabolic ID confinement) and (b) four different 
values of the width R^ of the channel (hard-wall confine- 
ment). 



tides diffusing in a channel of circular shape. In par- 
ticular, we will study the influence of the strength of the 
confinement (i.e., the depth of the potential profile across 
the channel) on the diffusive behavior. Without loss of 
generality, we will adhere to the specific conditions and 
parameters of the experimental set-up used in Ref . [25j . 
An additional advantage of this model is that the motion 
of the system of charged metallic balls [1^ is not over- 
damped, and we will solve the full Langevin equations of 
motion to study the diffusive behavior of the system. 

We consider N particles, interacting through a Yukawa 
potential [Eq. ([B])], which are embedded in a ring channel 
of radius rch ■ We define a parabolic confinement poten- 
tial across the channel in the form ([S]) where parameter 
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/? is chosen as follows: 



Vo — 



(a) 



'2Krrh sin 



in(^) 



2rch sin 



, (16) 



when all the particles are equidistantly distributed along 
the bottom of the circular channel. It should be noted 
that in this case, Vb is approximately equal to Vgs due 
to the weak Yukawa interaction, which slightly shifts the 
particles away from the bottom of the channel. Such a 
choice of Vq is related to the fact that we study the influ- 
ence of the confinement on the diffusion and, therefore, 
the potential energy of the particles must be of the or- 
der of the inter-particle interaction energy. Parameter 
To characterizes the distance where the external poten- 
tial reaches the value Vb/7, and Vgs is the energy of the 
ground state of the system of TV particles as defined by 
Eq. ([71) ■ Parameter 7 plays the role of a control param- 
eter. By changing 7 we can manipulate the strength of 
the confinement and, therefore, control the fulfillment of 
the single-file condition. Increase in 7 corresponds to a 
decrease in the depth of the confinement ([S]) which leads 
to the expansion of the area of radial localization of par- 
ticles. Therefore, an increase of 7 results in a similar 
effect (i.e., spatial delocalization of particles) as an in- 
crease of temperature, i.e., parameter 7 can be considered 
as an "effective temperature". Note that such a choice 
of the parameter that controls the confinement strength 
is rather realistic. In the experiment of Ref. [25| with 
metallic balls, the parabolic confinement was created by 
an external electric field, and the depth of the potential 
was controlled by tuning the strength of the field. 

To study diffusion of charged metallic balls, we solve 
the Langevin equation of motion in the general form (i.e., 
with the inertial term oc m). 



^~dt 



,) 



- VVcin) + F^, (17) 

where m = 2.5 x 10~^ kg [25| is the mass of a particle, 77 is 
the friction coefficient (inverse to the mobility). Here all 
the parameters of the system were chosen following the 
experiment and Xd = 4.8 x 10~^ m, F = 1 (which is 
a typical experimental value, see, e.g., also JlS]). Corre- 
spondingly, mass is measured in kg, length in m, and time 
in seconds. Also, following Ref. [25|, we took a channel 
of radius r^h = 9 mm (in the experiment [1^ , the exter- 
nal radius of the channel was 10 mm, and the channel 
width 2 mm; note that in our model we do not define the 
channel width: the motion of a particle in the transverse 
direction is only restricted by the parabolic confinement 
potential). We also took experimentally relevant num- 
ber of diffusing particles, N, varying from iV = 12 to 
= 40 (in the experiment [2^, the ring channel con- 
tained A = 12 or TV = 16 diffusing baUs). 

Fig. [7] shows the results of calculations of the trajec- 
tories of A = 20 particles diffusing in a ring of radius 
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FIG. 7: Trajectories of A'^ = 20 particles diffusing in a ring of 
radius rch = 9 mm for lO'' consequent time steps for different 
values of 7. 7=1 (a), 2 (b), 3 (c), 5 (d), 7 (e), 9 (f). 



Vch = 9 mm for the first 10^ MD steps for various values 
of the parameter 7. As can be seen from the presented 
snapshots, the radial localization of particles weakens 
with increasing 7. At a certain value of 7 this leads to 
the breakdown of the single-file behavior (Figs. [7l[c)-(f)). 



A. Breakdown of SFD 

It is convenient to introduce the distribution of the 
probability density of particles in the channel Pmd along 
the radial direction r. In order to calculate the function 
^rad('') we divided the circular channel in a number of 
coaxial thin rings. The ratio of the number of obser- 
vations of particles in a sector of radius to the total 
number of observations during the simulation is defined 
as the probability density Prad('"i)- Fig. [SI the proba- 
bility density Prad(?') is presented for different values of 
7. With increasing 7, the distribution of the probability 
density Prad('') broaden and the maximum of the func- 
tion Prad('') shifts away from the center of the channel 
(see Fig. [51). The latter is explained by the softening 
of the localization of particles with increasing 7, which 
tend to occupy an area with a larger radius due to the 
repulsive inter-particle interaction. Simultaneously, the 
distribution of the probability density Prad('') acquires 
an additional bump indicating the nucleation of a two- 
channel particle distribution . The observed broaden- 
ing and deformation of the function Prad('") is indicative 
of a gradual increase of the probability of mutual bypass 
of particles (i.e., the violation of the SF (single-file) con- 
dition, also called the "overtake probability" [43]) with 
increasing 7. 

Let us now discuss a qualitative criterion for the break- 
down of SFD, i.e., when the majority of particles leave 
the SFD mode. For this purpose, let us consider a parti- 
cle in the potential created by its close neighbor (which 
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FIG. 8: (Color online) The distribution of the probabil- 
ity density of particles Prad(7') in a circular channel of radius 
Vch = 9 mm along the radial direction r. The different curves 
correspond to various 7. Increasing 7 the width of the dis- 
tribution Pra.d{r) increases due to a weakening of the confine- 
ment. 




FIG. 9: (Color online) Spatial distribution of the poten- 
tial Vint (r, (j>) created by a particle (red (grey) circle) and the 
qualitative distribution of the probability density of particles 
in circular channel Prad(r) (green (light grey) line) along the 
radial direction r. The function Ar determines an approx- 
imate radial distance between particles when the potential 
barrier Utar becomes "permeable" for given temperature T. 
The function Araw characterizes a width of the distribution 
Pradir) at this temperature T. 



is justified in case of short-range Yukawa interparticle in- 
teraction and low density of particles in a channel) shown 
in Fig. m Different lines show the interparticle potential 
Vint as a function of angle </> for different radii r. For 
small values of 7, the center of the distribution Prad('') 
(see Fig. 7) almost coincides with the center of the chan- 
nel (i.e., with the minimum of the confinement potential 
profile) and the distribution Prad('') is narrow. There- 



fore, mutual passage of particles is impossible, i.e., the 
SF condition is fulfilled. The asymmetric broadening of 
the function Prad(^) with increasing 7 results in an in- 
creasing probability of mutual bypass of particles which 
have to overcome a barrier Ubar (see Fig. [9]). This be- 
comes possible when Ubar ^ ksT. In other words, the 
thermal energy ksT determines some minimal width Ar 
between adjacent particles when the breakdown of the 
SF condition becomes possible. 

It is clear that "massive" violation of the SF condition 
(i.e., when the majority of particles bypass each other) 
occurs when the halfwidth Arg^ of the distribution of 
the probability density Prad('') obeys the condition: 

Ars^u > Af . (18) 

The function Arg^, is defined by the ratio of the thermal 
energy k^T to the external potential Uconf{f) and is of 
the same order as Ar: 

• {ATl2f « keT. (19) 

Therefore the criterion (|18p can be presented in the form: 

Arsw « A^ > Af . (20) 

This qualitative analysis of the breakdown of the SFD 
regime clarifies the role of the width and the shape of the 
distribution of the probability density influenced by the 
asymmetry of the circular channel. 

B. Diffusion regimes 

The MSD {A(jP{t)) is calculated as a function of time 
t as: 

(A02(t)) = (j^^ J2 [^<^^^-(^ + ^) - A0„-(t)]M , 

\ JVpar-iVens - ■ j 

(21) 

where A'^par is the total number of particles of an ensem- 
ble and iVens is the total number of ensembles. In our 
calculations, the number of ensembles was chosen 100 for 
a system consisting of 20 particles. 

The time dependence of the MSD for different values 
of 7 is shown in Fig. [TUra)-(c). Initially the system ex- 
hibits normal diffusion, where (A0^) cx This regime 
is followed by an intermediate subdiffusive regime, where 
the (A0^) cx t" (0.5 < a < 1.0). For longer times, the 
system recovers "long-time" normal diffusion (see dis- 
cussions in Sec. IIIB), with (A0^) cx As in the 
case of straight channel geometry, this second crossover 
(i.e., from intermediate subdiffusion to "long-time" nor- 
mal diffusion) can also be due to two other reasons: (i) 
due to a collective (center-of-mass) diffusion or (ii) due to 
a single-particle jumping process. However, for the sim- 
ulations in the case of a circular geometry, the number 
of particles is relatively small (taking the fact that this 
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is a finite-size system), and therefore, the crossover from 
subhnear to linear regime is due to a collective (center-of- 
mass) diffusion. We further address this issue in Sec. V, 
where we consider a discrete site model and we exclude 
the center-of-mass motion. 

Fig. fTUTd) shows a as a function of 7. The function 
a{'j) experiences a monotonic gradual crossover from the 
a — 0.5 to a a < 1-regime. Note that the observed 
deviation from the normal diffusion behavior for large 
7 (Fig. [TUI is related to the presence of, though weak 
but nonzero, external confinement in the radial direction. 
This change of the diffusive behavior is explained by a 
weakening of the average radial localization of particles 
with increase of 7 (Fig. ^ and, as a consequence, by an 
increase of the probability of mutual bypass of particles. 

The observed crossover between the ID single-file and 
2D diffusive regimes, i.e., Q;(7)-dependence, shows a sig- 
nificant different qualitative behavior as compared to the 
case of a hard-wall confinement potential considered in 
Sec. lIIIi where a rather sharp transition between the two 
regimes was found [Fig. [IJd)]. The different behavior 
is due to the different confinement profiles and can be 
understood from the analysis of the distribution of the 
probability density of particles for these two cases. In the 
case of a hard- wall channel, the uncompensated (i.e., by 
the confinement) interparticle repulsion leads to a higher 
particle density near the boundaries rather than near the 
center of the channel (see Fig. [5] and Fig. IH^b)). As a 
consequence, the breakdown of the SF condition — with 
increasing width of the channel — happens simultane- 
ously for many particles in the vicinity of the boundary 
resulting in a sharp transition (see Fig. [Hd)). On the 
contrary, in the case of parabolic confinement, the den- 
sity distribution function has a maximum — sharp or 
broad, depending on the confinement strength — near 
the center of the channel (see Figs, ^a.) and|8]). With 
increasing the "width" of the channel (i.e., weakening 
its strength) , only a small fraction of particles undergoes 
the breakdown of the SF condition. This fraction grad- 
ually increases with decreasing strength of the confine- 
ment, therefore resulting in a smooth crossover between 
the two diffusion regimes. 




10 10" 10' 10' 10' 
t(s) 



DISCRETE SITE MODEL: 
LIMIT 



THE LONG-TIME 




The calculated MSD for different geometries and con- 
finement potentials allowed us to explain the evolution of 
the subdiffusive regime with varying width of the channel 
(or potential strength in case of a parabolic potential). 
However, the obtained results are only valid for the in- 
termediate regime and therefore they only describe the 
"onset" of the long-time behavior. The problem of ac- 
cessing the long-time behavior in a finite chain is related 
to the fact that sooner or later (i.e., depending on the 
chain length) the interacting system will evolve into a 
collective, or "single-particle" , diffusion mode which is 



FIG. 10: (Color online) (a)-(c): Log- log plot of the mean- 
squared displacement (MSD) {A0^) as a function of time for 
different values of the "effective" temperature 7 = (a) 1, (b) 
2, and (c) 3. Here {N^ns = 100, A^par = 20). (d) The diffu- 
sion exponent q as a function of 7. Increase of the "effec- 
tive" temperature 7 leads to the gradual transformation of 
the single-file regime of diffusion into the diffusion regime of 
free particles. 



11 



characterized by a = 1.0. Thus the question is whether 
the observed behavior holds for the long-time limit, i.e., 
is the transition from to t^-^ behavior smooth? 

To answer this question, we considered a simple model, 
i.e., a linear discrete chain of fixed sites filled with either 
particles or "holes" (i.e^sites not occupied by particles) 
(for details, see Ref. (Sg); this model was also recently 
used in Ref. [13] ) . The particles can move along the chain 
only due to the exchange with adjacent vacancies (i.e., 
with holes). Within this model, the long-time diffusion 
behavior was described analytically for an infinite linear 
chain as well as for a finite cyclic chain '38^. In partic- 
ular, this model predicts that: (i) If the chain is infinite 
then the long-time power law of the diffusion curve a is 
0.5 (i.e., MSD (Aa;2(i)) « t°-^); (ii) If the chain is finite 
then the subdiffusive regime with a = 0.5 is followed by 
either a = 1.0 regime (if the cyclic boundary condition 
is realized), or by a = regime, i.e., the regime of satu- 
ration (if no cyclic boundary condition is imposed [l^). 
The latter regime is reached for times longer than the 
"diffusion time" of a "hole" along the whole chain ichain- 

Let us now apply this model to a finite-size chain of 
particles. For this purpose, we assume that adjacent par- 
ticles are able to exchange their positions with some prob- 
ability P at every time step. For example, probability 
P = 0.1 means that a couple of any adjacent particles 
certainly exchange their positions once for every 10 time 
steps. 

The results of our calculations of the MSD performed 
using this model are presented in Fig. lTTT a). We used the 
following parameters: the chain length is Ng = 150 sites 
and Nh = 1 hole. Averaging was done over 1000 ensem- 
bles. The calculation was performed for the following val- 
ues of the probability: P = 0, 10~^ 10"'', lO^^, Q.Ol, 0.1, 
and 1. 

We see in Fig. [TlT a') clearly the above-mentioned two 
diffusion regimes, i.e., with the MSD {/S.x^{t)) oc t^-^ 
and oc i^". The characteristic time ichain shifts towards 
lower values with increasing P. However this analysis 
(Fig. [TlTa')') does not allow to distinguish the contri- 
butions to the long-time behavior (oc t^ °) due to: (i) 
the breakdown of single-file condition (i.e., diffusion due 
to particle exchanges), and (ii) the "collective" diffusion 
(chain "rotation"). 

To overcome this difficulty, we exclude the "collective" 
diffusion of the system and introduce a modified MSD 
(Aa:^(t))corr (which is so-called "r oug hness" of the system 
of particles, as discussed in Ref. |28'|) as follows: 

(Aa;2),„„ = {{x-xf), 

where (...) is the average over time; x is the average of an 
ensemble of particles at a given time, or "collective" co- 
ordinate. It should be noted that {x) ^ x. If the system 
does not experience "collective" diffusion then x{t) = 
and the modified MSD coincides with the conventional 
one: 

(Ax2),o„ = (a;2). 




t (steps) 




t (steps) 

FIG. 11: (Color online) Log-log plot of the MSD (Aa;^(t)) (a) 
and corrected MSD (Aa::^(f))corr (b) as a function of time for 
different values of the probability P of bypassing. Averaging 
was done over A'sim = 1000 ensembles. 



The diffusion curves calculated by using the modified 
MSD are presented in Fig. HHb). For P = 0, the diffu- 
sion curve (shown by black open squares) after the sub- 
diffusive regime reaches saturation (i.e., (Aa;^(t))corr = 
const). The observed behavior is similar to that of a 
finite linear chain with fixed ends (see Ref. [l3|)- For 
P 7^ 0, all the diffusion curves in the long-time limit are 
characterized by a = 1.0, independent of the value of the 
probability P, as seen in Fig. inT b). In other words, the 
long-time diffusion does not depends on the probability 
of mutual exchanges of particles and has the same long- 
time behavior for any probability P 7^ 0. Here we would 
like to emphasize again that the long-time behavior of 
the diffusion curves is free from the "collective" diffusion 
effect and is only determined by particle jump diffusion. 
Increasing a number of sites in the model corresponds, 
in fact, approaching to the model of infinite chain. We 
have found that the increasing a number of sites leads 
to growth of the (Ax^(i))corr limit of saturation, on the 
one hand, and to a shift of tchain to larger on the other 
hand. Hence, extrapolating our results to the case of infi- 
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nite chain, we can conclude that in this case as well as in 
the case of finite-size chain, the breakdown of single-file 
condition leads to an abrupt transition from subdiffusive 
to the normal diffusion regime. 

The difference in the diffusive curves is just the time 
Ttran from subdiffusivc regime to the normal regime: for 
low P it (rtran) IS loug cuough while for high P it (rtran) 
is short. It is easy to see that rtran ~ 1/P{%)- Thus, we 
can conclude that in the long-time limit the transition 
from t'^-^ to ° behavior is abrupt. Note that our cal- 
culations performed using the modified MSD (Ax^(t))corr 
reproduce the results of Ref. jl^] for a closed "box" . This 
is explained by the fact that in the closed "box" geom- 
etry the center of mass (or collective) diffusion is zero, 
and it is natural that the roughness (see Ref. [1^) and 
the particles diffusion coincide. 

VI. CONCLUSIONS 

We have studied a monodisperse system of interact- 
ing particles subject to three types of confinement po- 
tentials: (i) a ID hardwall potential, (ii) a ID parabolic 
confinement potential which both characterize a quasi- 
ID system, and (iii) a circular confining potential, which 
models a finite size system. In order to study the diffusive 
properties of the system, we have calculated the mean- 
squared displacement (MSD) numerically through molec- 
ular dynamics (MD) simulations. For the case where 
particles diffuse in a straight line in a QID channel, dif- 
ferent diffusion regimes were found for different values 
of the parameters of the confining potential {x or Rw)- 
We have found that the normal diffusion is suppressed 
if the channel width is between 0.20 and 0.50 (or by 
2.0 < X < 3.5, for the case of parabolic ID confinement), 
leading the system to a SFD regime for intermediate time 
scales. For values of Rw ^ 0.56, particles will be able to 
cross each other and the SFD regime will be no longer 
present. 

The case of a circular channel corresponds to, e.g., the 



set-up used in experiments with sub-millimetric metallic 
massive balls diffusing in a ring with a parabolic potential 
profile created by an external electric field. The strength 
of the potential (which determines the effective "width" 
of the channel) can be tuned by the field strength. Con- 
trary to the case of hard- wall confinement, where the 
transition (regarding the calculation of the scaling expo- 
nent (a) of the MSD (Ax^(t))oc t") is sharp, a smooth 
crossover between the ID single-file and the 2D diffusive 
regimes was observed. This behavior is explained by dif- 
ferent profiles for the distribution of the particle density 
for the hard-wall and parabolic confinement profiles. In 
the former case, the particle density reaches its maximum 
near the boundaries of the channel resulting in a massive 
breakdown of the SF condition and thus in a sharp transi- 
tion between the different diffusive regimes. In the latter 
case, on the contrary, the density distribution function 
has a maximum near the center which broadens with de- 
creasing strength of the confinement. This results in a 
smooth crossover between the two diffusion regimes, i.e., 
SFD and 2D regime. The analysis of the crossing events, 
i.e., the rate of the crossing events Wc as a function of 
the confinement parameter % or R^,, supports these re- 
sults: the function uJc{x/Rw) displays a clear signature 
of either "smooth" or "abrupt" behavior. 

We also addressed the case of a finite discrete chain 
of diffusing particles. It was shown that in this case 
the breakdown of the single- file condition (i.e., when the 
probability P of particles bypassing each other is non- 
zero) leads to an abrupt transition from a subdiffusive 
regime to the normal diffusion regime. 
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